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Abstract 

Since its decovery during the late 90's, the dimming of distant SN la 
apparent luminosity has been mostly ascribed to the influence of a mys- 
terious dark energy component. Based upon the cosmological "principle" 
hypothesis, this interpretation has given rise to the "concordance" model, 
developed in the context of a Friedmann-Lemaitre cosmology. However, a 
caveat of this reasoning is that the cosmological "principle" derives from 
a philosophical Copernican assumption and has never been tested. Fur- 
thermore, a weakness of its conclusion, i. e., the existence of a negative- 
pressure fluid or a cosmological constant, is that it would have profound 
implications for the current theories of physics. This is why we have 
proposed a more conservative explanation, ascribing the departure of the 
observed universe from an Einstein-de Sitter model to the influence of 
inhomogeneities. This idea has been independently developed by other 
authors and further enlarged to the reproduction of different cosmological 
data. We review here the main proposals which has been made along 
these lines of though and present some prospects for future developments. 

1 Introduction 



Since its decovery during the late 90 's, the dimming of distant SN la apparent 
luminosity has been mostly ascribed to the influence of a mysterious dark en- 
ergy component PQ which is supposed to drive an acceleration of the Universe 
expansion. 

Based upon the cosmological "principle" hypothesis and coupled to an anal- 
ysis of the other cosmological data, this interpretation has given rise to the "con- 
cordance" or ACDM model, developed in the context of a Friedmann-Lemaitre 
cosmology. A drawback to this reasoning is that the cosmological "principle" 
merely derives from a philosophical Copernican assumption which still deserve 
to be tested 0[3]. 
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The dark energy component, a cosmological constant or a negative pressure 
fluid, represents, in the "concordance" model, about 73% of the Universe energy 
density. However, a cosmological constant is usually interpreted as a vacuum 
energy of which particle physics cannot explain such an amplitude and a neg- 
ative pressure fluid remains a mysterious phenomcnum. This is known as the 
cosmological constant problem. 

Another feature of the supernova data analysed in a Friedmannian frame- 
work is to yield a late-time acceleration of the expansion rate, about the epoch 
when structure formation enters the nonlinear regime. This would imply that 
we live at a time when matter density energy and dark energy are of the same 
order of magnitude. This is known as the coincidence problem. 

General relativity having only been tested up to scales of order a plane- 
tary system, a second type of explanation has been proposed which implies a 
modification of this theory at larger distance scales [I]. 

Our current purpose is to review the works dedicated to a more simple 
and natural proposal, which makes only use of known physics and phenomena. 
Since it appears that the onset of the apparent acceleration and the beginning 
of structure formation in the Universe are concomitant, the idea that the SN 
la observations could be reproduced by the effect of inhomogeneities have been 
put forward. This interpretation has been first proposed independently by a few 
authors [21 13 [7] , shortly after the release of the data. Then, after a period 
of relative disaffection, it has experienced a reniewed interest, specially during 
the last two years. 

Another natural explanation is the possible effect of an actual geometrical 
cosmological constant. Its value such as it appears in the ACDM model has been 
indeed predicted in the framework of scale relativity, as A = 1.36 10~ 56 cm~ 2 , 
long before the release of the supernova results [H [9]. However, we shall not 
develop this approach here since it is beyond the scope of this contribution. 

2 Presentation and discussion of the methods 

It is not obvious that one can keep using FLRW models to interpret the high 
precision observational data collected in a regime where most of the mass is 
clumped into or is forming structures. Three main physical effects might actu- 
ally be missing [10 : (i) the overall (average) dynamics of such a universe could 
be significantly different from the FLRW one. This effect is usually taken into 
account by calculating so-called backreaction terms added to effective equations 
for the dynamical evolution of the physical quantities under study, (ii) light 
propagation in a clumpy universe might be different from that in a homoge- 
neous one and the luminosity distance-redshift (LDR) relation could hence be 
affected, (iii) the fact that we have only one observer could influence the re- 
sults, since we might live in an under or over-dense region which should induce 
significant corrections in the data interpretation. 

What is the most subtle point of this issue is that what is observed in the 
supernova data is not an accelerated universe expansion (this is only an artefact 
due to the a priori assumption that our universe, even in the local region of 
structure formation, can be represented by a FLRW model) but only a dimming 
of the SN la luminosity taken into account by the infered LDR relation. 

However, many authors have focussed their attention on effect (i) and back- 
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reaction has been studied with two types of methods: (a) in the linear regime of 
small amplitude fluctuations, perturbative expansions which were subsequently 
averaged out, (b) for the more general problem of the dynamics induced by all 
scale inhomogeneities including the nonlinear regime, spatial averaging of non- 
perturbed models. The local dynamics has also been analysed using different 
exact solutions of the general relativity equations. 

Effects (ii) and (iii) have been mostly studied using toy models constructed 
with such a kind of exact nonaveraged inhomogeneous solutions. 

2.1 Perturbative analysis 

The perturbative analysis is a method employed when the deviations from ho- 
mogeneity and isotropy are assumed to be "small" (see Mukhanov et al. [11] 
for a review). 

One considers two spacetimes, the physical, perturbed spacetime and a fic- 
tious background spacetime described by a FLRW model. A one-to-one corre- 
spondence between points in the background and points in the physical space- 
time carries the coordinates labelling the points in the background over into 
the physical spacetime and defines a choice of gauge. The perturbation in some 
quantity is the difference betwen the values it possesses at a point in the physical 
spacetime and at the corresponding point in the background spacetime. Since 
we arc interested in considering the influence of inhomogeneities, the perturbed 
quantity is here the density. 

To study the influence of such fluctuations on the expansion rate of the Uni- 
verse, one identifies a local physical variable which describes this expansion rate, 
calculates the backreaction of the cosmological perturbations on this variable 
to a given order, and then spatially averages the result. Note that it is of the 
upmost importance to avoid the deficient procedure consisting of calculating 
an "observable" from the spatially averaged metric, which, in general, does not 
give the same result as calculating the spatial average of the observable [12] [T3] . 
As we shall stress in Sec. 12.21 it is also crucial to define the hypersurfaces on 
which the averaging is performed by a clear physical prescription. 

The cleanest available calculation of the effect of density fluctuations on the 
averaged expansion rate of a matter-dominated universe up to second order in 
the metric variables has been performed by Kolb et al. [14j for an application 
to superhorizon scale perturbations, in the adiabatic case. Owing to the per- 
turbation development employed, this method is only consistent for the study 
of the effect of very large scale fluctuations, with small amplitude, e. g., in 
the linear regime of structure formation [HI [H] \W\ . It is also sensitive to the 
order at which the development is performed. It has been demonstrated by 
Notari [17] that, while at early times the contribution of subhorizon inhomo- 
geneous gravitational fields is perturbatively subdominant, the series is likely 
to diverge around the time of structure formation. There one must employ a 
nonperturbative method. 

To deal with this issue either spatial averaging of a nonperturbed solution 
or exact solutions of Einstein's equations are usually used. A new formalism 
for the study of nonlinear perturbations in cosmology, based on a covariant and 
fully nonperturbative approach, has recently been proposed [18j but it has not 
yet been applied to our purpose. 
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2.2 Spatial averaging of nonperturbed models 

Spatial averaging aims at obtaining the impact of a given inhomogeneity profile 
upon the assumed large scale FLRW background, in terms of backreaction terms 
added to the Friedmannian evolution equations of scalar quantities such as the 
expansion rate, the energy density, the isotropic pressure. These backreaction 
terms can be assimilated to a black energy component if they exhibit the right 
properties. 

In general relativity, spatial averaging is very much involved |19j since the 
equations which determine the metric tensor and the quantities calculated from 
it are highly nonlinear. However, when modelling the Universe, the usual 
method is to use continuous functions representing, e. g., energy-density, pres- 
sure, or other kinematical scalars of the velocity field, implicitly assuming that 
they represent volume averages of the corresponding fine-scale quantities. But 
we know that our local Universe is highly inhomogeneous from the scales of 
planetary systems up to currently unknown distances 20J. Anyhow, the scale, 
i. e., the size of the volume, over which the averagings are performed are never 
explicitly defined, while the results of this process obviously depend on it. 

Moreover, a volume average is a well-defined quantity for scalars only. For 
vectors, and all the more for tensors, it leads usually to noncovariant quanti- 
ties. Another drawback is that a gauge problem arises when relating the "true" 
and the averaged metric. Scalar quantities only are invariant under coordinate 
transformations, not tensors. One must also be aware that, in a generic space- 
time, there are no preferred time-slice one could average over and the results 
depend on the choice of the hypersurfaces on which this average is performed 

But the main issue is the noncommutating property of the two operations: 
averaging the metric and calculating the Einstein tensor. In other words, 
the Einstein tensor calculated from an averaged metric and energy-momentum 
tensor is not equal to the Einstein tensor first calculated from the fine-scale 
metric and energy-momentum tensor, then averaged. However, in the stan- 
dard cosmological approach, the Universe is modelled by adopting exactly the 
wrong method: take a metric which is assumed already averaged, calculate 
the corresponding Einstein tensor and equate it to an already averaged energy- 
momentum tensor. 

The question raised by the method consisting of determining the parameters 
of an a priori assumed FLRW model from observational data is the "fitting prob- 
lem" [22] . The departure of the "real" Universe from the averaged one is known 
as the "backreaction" effect. The study of backreaction can be implemented 
by two methods: either one tries to obtain directly the equations satisfied by 
the averaged quantities, with minimum assumptions as regards an underlying 
background, or one takes as a background a FLRW model and analyses the 
effect of linear perturbations applied on it. In this section, we consider the first 
approach. The second one has already been studied in Sec. 12.11 

This issue has been delt with by Buchert and collaborators for simplified 
inhomogeneous cosmological models with an irrotational perfect fluid as the 
gravitational source. First, Buchert and Ehlers [23| have proposed a simple 
averaging procedure which have been used by Palle [24] to deal with the cos- 
mological constant problem. Then, Buchert [551 US] has developed another 
procedure aimed at constructing an "effective dynamics" of spatial portions of 
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the Universe from which observable average characteristics can be infered like 
the Hubble constant, the effective 3-Ricci scalar curvature and the mean density 
(and isotropic pressure) of a spatial domain. 

Relations between average scalar sources (energy density, pressure) and an 
average scalar geometric quantity, the expansion rate, have been derived for an 
adapted foliation of spacetime. They involve domain dependent backreaction 
terms that have been splitted into a "kinematical" backreaction comparing the 
variance of the expansion rate to the shear, and a "dynamical" backreaction, i. 
e., pressure forces. These equations show that the averaged shear fluctuations 
tend to increase the expansion rate as do the averaged energy source terms 
(provided the energy condition holds) , while the averaged expansion fluctuations 
have an opposite effect and therefore work to a stabilization of structures. The 
dynamical backreaction can do both. 

However, an averaging procedure is not complete unless one also averages the 
geometrical inhomogeneities. Since geometrical fields are tensorial variables for 
which possible strategies of averaging are not straightforward, Buchert and Car- 
fora [57] have suggested a Lagrangian smoothing of these variables as opposed 
to their Eulerian averaging on spatial domains. Curvature fluctuations turn out 
to be crucial and may even outperform the effect of kinematical fluctuations. 
The authors define effective cosmological parameters that would be assigned 
to the smoothed cosmological spacetime. These parameters are "dressed" after 
smoothing out the geometrical fluctuations. Relations between the "dressed" 
and "bare" parameters are derived. The former provide the framework for inter- 
preting observations with a "Friedmann bias", i. e., as if the observer was living 
in a Friedmannian universe. The latter represent the actual inhomogeneous 
cosmological model, spatially averaged. 

In subsequent articles [28j [29] , the same authors have identified two effects 
that quantify the difference between "bare" and "dressed" parameters: the "cur- 
vature backreaction" and the "volume effect". To summarize we can say that, 
in the smoothed model, the averaged scalar curvature is dressed both by the 
volume effect and the curvature backreaction effect. The volume effect is ex- 
pected when comparing two regions of distinct volumes, but with the same 
matter content, in a constant curvature space. The backreaction term encodes 
the deviation of the averaged scalar curvature from a constant curvature model, 
e. g., a FLRW space section. 

However, the interpretation of cosmological parameters remains far from 
trivial. One must also average on the observer's light cone in which case the 
above effects interact with the time evolution of the model. Moreover, smooth- 
ing will have to be performed in a dynamical setting. Apostolopoulos et al. [30] 
have uncovered another subtelty of inhomogeneous cosmology: the volume in- 
crease at a given point results from averaging over various directions a possible 
anisotropic acceleration parameter and might thus not be the most appropriate 
for an expansion characterization. 

2.3 Use of exact solutions of Einstein's equations 

The method consisting in using exact solutions to modelize the inhomogeneities 
observed in the Universe is the most straightforward and devoided of theoretical 
pitfalls. It is adapted to represent either strong or weak inhomogeneities. 
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For mathematical simplification and also to account for the local quasi- 
isotropy of the CMB as measured on our wordline, most of the retained models 
exhibit spatial spherical symmetry. Therefore, some authors have claimed that 
unphysical properties of this symmetry, and specially of the Lemaitre-Tolman- 
Bondi (LTB) class [31] [32], might prevent these models to be used in this frame- 
work, specially when the observer is assumed to be located at the centre. How- 
ever, in the literature, the observer has been put either at the center of the 
model j33[ [33] or offcenter [35]. Moreover, a spherically symmetric model can 
represent an anisotropic inhomogeneous universe averaged over angular scales, 
which is not physically worse than a uniformly averaged universe such as in the 
FLRW picture. 

Other spherically symmetrical models has been used in this framework, i. 
e., peculiar classes of Stephani solutions. One feature of the Stephani models 
that has been the subject of much debate is their matter content [31? ] \W \ 137]. 
Individual fluid elements can behave in a rather exotic manner, e. g., exhibiting 
negative pressure. But this is also the case of dark energy. 

Avoiding spherical symmetry, Moffat [35] has considered a peculiar class of 
models of the Szafron family. Other exact solutions have consisted of FLRW 
patch(es) embedded in a FLRW background with different energy densities 

Bans. 

Depending on the authors, these exact solutions have been used either to 
barely fit the data as they are obtained by observation or to provide the in- 
homogenous models on which an averaging procedure "a la Buchcrt" is imple- 
mented [IT] . 

3 Studied physical quantities 

3.1 The deceleration parameter 

When reasoning in the framework of a Friedmannian cosmology, the dimming 
of the supernovae is associated with an acceleration of the Universe expansion. 
This is why a number of authors have focussed on the issue of either demon- 
strating or ruling out an effect of the inhomogeneities on the expansion rate. 

Some have tried to derive [HI [33J or rule out [3H1 S3] no-go theorems. How- 
ever, when spatially averaged, a physical quantity associated with the expansion 
rate behaves quite differently [331 [33 HB] ■ It is thus difficult to yield general 
rules from such theorems. 

Some have stressed that the definition of a deceleration parameter in an in- 
homogeneous framework is tricky. Hirata and Seljak [42! have proposed four 
different definitions of such a parameter. Apostolopoulos et al. [3U] have shown 
that an observer located away from the centre of a spherically symmetric config- 
uration can measure acceleration or deceleration along radial or perpendicular 
directions depending on local under or over-density. This demonstrates that, 
even locally, the effect of inhomogeneities on the dynamics of the Universe is 
not trivial. 

Ishibashi and Wald [37J have also argued that an averaged quantity repre- 
senting the scale factor or the deceleration parameter may accelerate without 
there being any observable consequence. 

Another unexpected effect has been put forward by Tomita [7j [39l [40] who 
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has considered a cosmological model composed of a low-density inner homoge- 
neous region connected to an outer homogeneous region of higher-density. Both 
regions decelerate, but, since the void expands faster than the outer region, an 
apparent acceleration is experienced by the observer located inside the void. 

We can therefore conclude that the computation of some local quantity (gen- 
erally the deceleration parameter), eventually subsequently averaged, and be- 
having the same way as in FLRW models with dark energy can lead to spurious 
results l.'i ;48J [3H| and must therefore be avoided. 

3.2 The Effective Stress Energy Tensor (SET) 

Martineau and Brandcnberger [SU] have tried to estimate the effect of a back- 
reaction of Super-Hubble modes, by computing it in terms of the second order 
Effective SET of cosmological perturbation theory. This has been criticized by 
Ishibashi and Wald 47] who have argued that a large SET implies the contribu- 
tion of higher order terms in the perturbative scheme, therefore impairing the 
claimed results. 

3.3 The LDR relation 

The LDR relation is the only direct product of the supernova data, obtained 
without any a priori cosmological assumption. This relation is therefore the 
best observable to be fitted and this has been the aim of a number of works 

0E1IIDII3DII3DI331IMI2ZII5I1I521I3I531. 
4 Main physical results 

It has been proposed that energy density perturbations of wave-length larger 
than the Hubble radius, generated during inflation, might produce upon cosmic 
parameters effects that could mimic an accelerated expansion of the Universe 
[141 I54[ I15[ 150] . But it has subsequently been shown by a number of authors 
[42], 143], 155] 156] [57] , among whom some of the proponents themselves [58] [59] , 
that this could not be the case. 

Working in the framework of a peculiar toy model, Biswas et al. [110 have 
computed a minimal overall (average) effect which amounts to a correction in 
apparent magnitudes at all rcdshifts of order Am ~ 0.15. 

This effect is to be added to the estimation made by Buchert and Carfora 
[29] of the volume effect alone in a naive swiss-cheese model. Its magnitude 
of 67% reduces the necessary dark energy in the concordance model roughly 
from 70% to 50% [BU] . A mismatch of similar magnitude has been reported by 
Hellaby |JT], using volume matching [22] in LTB models of clusters and voids. 

Studies of the constraints which apply to a dust model of universe if one 
wants to explain dark energy by the backreaction effect of inhomogeneities have 
shown that a negative average spatial curvature is mandatory to compensate a 
strong backreaction increasing the expansion rate of the model [56l[6TH[2Tl[S2l[551 
[3H @S] . This result is actually consistent with the fact that most of the exact in- 
homogeneous models which reproduce best the observed LDR relation are those 
which exhibit negative spatial curvature near the observer, or, equivalcntly, 
where the observer is located in an underdense region [TOj, 134] [39] [52j [53] [64] . 
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The assumption that we might be located within an underdense region seems 
to be consistent with observations leading to the identification of a Local Void 
and of its suggested expansion [301 [SS] . 

The matching of the observed LDR relation has been successfully performed 
by a number of proposed models: homogeneous void models [39j . LTB models 
with a centered observer [3U [HS [52] , LTB models with an outcentre observer 
[TUl IM] , Stephani models [37J HE EH 133] • Most of these models reproduce the 
ACDM LDR relation up to a redshift of the order unity. This feature, which 
has been considered as a ruling out drawback by some authors [HU 132] , must be 
viewed as a nice way out of the coincidence problem. The physical explanation 
is that the appearance of small scale inhomogeneities corresponds to the onset of 
structure formation at the time (around z ~ 1) where apparent "acceleration" 
begins to be observed |3T] . 

However, it has been shown [3J that the problem of fitting a given LDR 
relation with a peculiar inhomogeneous (LTB) model is completly degenerate. 
Therefore, the model parameters must be constraint by other cosmological data. 
This has been done by some authors: the model proposed by Biswas et al.fTU] 
reproduces matter abundance, the first acoustic peak in the CMB power spec- 
trum and the baryon oscillations; Alnes et al. [31] obtain the matter density 
measured at low z and the two first CMB peaks; Godlowski et al. [33], the 
three first CMB peaks; Bolejko [31] uses observed cosmological properties to 
constraint his models. The main challenge remains to fit the whole set of avail- 
able cosmological data with a given model. Since very few exact solutions to 
Einstein's equations can be of use in a cosmological framework, only oversimple 
toy models have been studied up to now to deal with this issue. However, a 
project to begin implementing this is currently underway |67j . 

5 Conclusion and prospects 

One can find in the literature inhomogeneous cosmological toy models able to 
solve both the cosmological constant and coincidence problems, i. e., to mimic 
an "accelerated expansion" with no need for dark energy up to the epoch when 
structure formation enters the nonlinear regime, i. e., around z ~ 1. 

It has been shown that inhomogeneities likely to solve these problems must 
be of the subhorizon and strong type, which cannot be studied with perturbation 
methods. Averaging and smoothing procedures have been proposed which can 
provide some insight into the issue for very simple cases but which much be 
used with care since they are not devoided of pitfalls and incompleteness. 

Exact solutions of Einstein's equations have the nice property of being able 
to modelize both strong and weak inhomogeneities. What these models must 
reproduce is not an accelerated expansion, which is an artefact of the homoge- 
neous assumption, but the observed dimming of the SN la luminosity, i. c., 
the LDR relation. Some classes have been nicely fitted to this relation at low 
redshifts. However, the proposed toy models do not pretend to be fair represen- 
tations of our neighbouring patch of universe. Actually, some of them exhibit 
controversed properties issued from their spherical symmetry (LTB models) or 
their matter content (Stephani models). But since very few exact solutions to 
Einstein's equations can be of use in a cosmological framework, these simple 
models are usefull to obtain some insight into the issue. 
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Even if some of them have been shown to match other cosmological obser- 
vational constraints, the main challenge remains to fit the whole set of available 
cosmological data. A project to begin implementing this issue is currently un- 
derway [H7] . A solution, if any, would be to use non pathological exact inhomoge- 
neous solutions, reproducing the nearby Universe, coupled to or asymptotically 
matching nearly homogeneous ones valid up to last scattering. 
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